A transparent expression of the ^-Condensate's 

renormalisation 



Ph. Boucaud Q , F. De Soto b , A. Le Yaouanc a , J. P. Leroy a , J. Micheli a , H. Moutarde c , O. Pene a , 

J. Rodriguez Quintero rf 

a Laboratoire de Physique Theorique * 
Universite de Paris XI, Batiment 210, 91405 Orsay Cedex, France 
b Dpto. de Fisica Atomica, Molecular y Nuclear 
Universidad de Sevilla, Apdo. 1065, 41080 Sevilla, Spain 
c Centre de Physique Theorique Ecole Polytechnique, 91128 Palaiseau Cedex, France 
d Dpto. de Fisica Aplicada e Ingenieria electrica 
E.P.S. La Rabida, Universidad de Huelva, 21819 Palos de la fra., Spain 

We give a more transparent understanding of the vacuum expectation value of the renormalised 
local operator A by relating it to the gluon propagator integrated over the momentum. The 
quadratically divergent perturbative contribution is subtracted and the remainder, dominantly due 
to the 0(l/p 2 ) correction to the perturbative propagator at large p 2 is logarithmically divergent. 
This provides a transparent derivation of the fact that this 0(l/p 2 ) term is related to the vacuum 
expectation value of the local A 2 operator and confirms a previous claim based on the operator 
product expansion (OPE) of the gluon propagator. At leading logarithms the agreement is quanti- 
tative, with a standard running factor, between the local A 2 condensate renormalised as described 
above and the one renormalised in the OPE context. This result supports the claim that the BRST 
invariant Landau-gauge A 2 condensate might play an important role in describing the QCD vacuum. 
P.A.C.S.: 12.38.Aw; 12.38.Gc; 12.38. Cy; 11.15.H 

I. INTRODUCTION 

In a series of lattice studies [1-4] the gluon propagator in QCD has been computed at large momenta, and it was 
shown that its behavior was compatible with the perturbative expectation provided a rather large 1 /p 2 correction was 
considered. In an OPE approach this correction has been shown [2,3] to stem from an A 2 gluon condensate which 
does not vanish since the calculations are performed in the Landau gauge. It was also claimed [4] that this condensate 
might be related to instantons. 

The role of such a condensate in the non-perturbative properties of QCD, in particular its relation to confinement, 
has been studied by several groups [5,6]. Of course any physics discussion about the A 2 condensate necessitates a clear 
definition of what we speak about, i.e. it needs a well defined renormalisation procedure to define the renormalised 
local A 2 operator, since j4(0) 2 is a quadratically divergent quantity as can easily be seen in perturbation theory. 
A renormalisation of A 2 was defined in [2,3] within the OPE context which we now briefly summarise. It uses the 
notion of "normal order product" in a "perturbative vacuum" which is annihilated by the fields A. It implies that 
<: A(0) 2 :>p 0r t= in the perturbative vacuum^. The contribution to <: A(0) 2 :> in the true QCD vacuum is then 
of non-perturbative origin. It has only logarithmic divergences and it is multiplicatively renormalised. Of course this 
notion of a perturbative vacuum in which Fock expansion could be performed has not a very transparent physical 
meaning especially in a non-perturbative context such as the numerical Euclidean path integral method. 



* Unite Mixte de Recherche du CNRS - UMR 8627 

^The symbol ": ... :" represents the normal ordered product in this perturbative vacuum. 



1 



A 2 is not a gauge invariant operator but the bare A 2 condensate is a very special object since, by definition, it is a 
minimum of the gauge orbit [6]. In other words, some important physics seems to lie beneath the BRST invariance 
of A 2 in Landau gauge. The authors of ref. [5] discussed on the generalised* composite operator A^A^ + 2iacc, 
which is BRST invariant in the manifestly Lorentz covariant gauge, and examined the survival of this invariance after 
renormalisation. In this note, although in a different context, we also examine the same point: the subtle relationship 
between the minimum of bare A 2 in the gauge orbit and any gauge-independant physical phenomenology associated 
to the renormalised condensate [6], emerging for instance from the OPE analysis [1-4]. To this aim, we will derive the 
renormalised A 2 vacuum expectation value without using the normal ordering but using only the OPE expansion of 
the gluon propagator § . It will provide a more transparent definition, related directly to a quantity which is actually 
measured. 

We start from the observation that the non renormalised < A(0) 2 > is related to the integral of the gluon propagator 
over momentum. Hence it is expected that the non-perturbative contribution to A 2 has to do with the non-perturbative 
contribution to the gluon propagator. The latter contains precisely 1/p 2 contributions due to the A 2 condensate at 
large momenta, and also strong deviations from perturbative QCD at small momenta, see fig. 1 (taken from [1-4]). 
How does this fit together? 



II. BARE, PERTURBATIVE AND NON-PERTURBATIVE A 2 



It is possible in principle from lattice calculations to define the non-perturbative gluon propagator in the Landau 
gauge. Lattice calculations provide the bare gluon propagator. From the gluon propagator computed with a series 
of different values of the lattice spacing one can in principle compute the renormalised gluon propagator from zero 
momentum up to as large a momentum as one wishes. An example of such a non-perturbative propagator is shown 
in fig. 1. We can choose for example the MOM renormalisation scheme**, such that 



GkV = M a ) = 



1 



(1) 



This implies a renormalisation of the gluon fields 



A v r = Z 3 (n) iA v bare , Z 3 (n) = /j 2 G 



(2) 
bare 



(2) 



The renormalisation constant Z% has to be understood as related to any regularisation method and any value of the 
UV regulator provided that the latter is larger than the momenta carried by the gluons. The coupling constant is 
also renormalised in the MOM scheme. The Yang-Mills theory is thus fully renormalised and from now on we will 
consider only renormalised gauge fields and propagators. 
The propagator is defined in Euclidean space by 



J < A- K (0)Ai R (x) >=5 afi 

Inverting the Fourier transform, 



, _ Pj^Pu 
W p 2 



Gr( P 2 ) 



E < ^r(oms r (o) >= / rfV4 2 V) = 3 -^A J P V<4V) 



(3) 



(4) 



* Landau gauge is recovered in the limit a — > 

§ Of course the normal ordering has been used in [2,3] to compute the anomalous dimension of A 2 and the Wilson coefficient 

C2 

** Notice that the chosen renormalisation scheme is not relevant in our argument in this paper, but we clearly need a scheme 
in which non-perturbative quantities coming from lattice simulations can be accommodated. MOM is one of the simplest. On 
the contrary the MS scheme does not fulfill this condition. 
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This integral is quadratically divergent in the ultraviolet. Indeed, if the gauge fields and the coupling constant have 
been renormalised, the local A 2 operator has not yet. Let us introduce an ultraviolet cut-off A and define 

< (AMf >a= KN l~, 1] £ P 2 dp 2 G%\pi). (5) 

where (Ar(/j,)) 2 refers to the square of the gauge fields renormalised at the scale fi, but where A 2 has not been 
renormalised as a local product of operators. The symbol "< ... >" represents the vacuum expectation value. (Ar) 
is clearly an UV divergent quantity. The index A refers to the ultraviolet cut-off and /x to the renormalisation point 
for the gauge fields and the coupling constant. The cut-off A has nothing to do with the lattice cut-off a -1 . The 
renormalisation in eqs. (1) and (2) has eliminated any dependence in the different lattice spacings which have been used 
to produce the renormalised propagator. A is introduced simply to control the quadratic and logarithmic divergences 
we encounter here. 

The dominant contribution to this integral is the perturbative one. To separate the perturbative contribution from 
the non perturbative we will now use the results of [3] 

i r (2) ( 2 2,_ P 2 G (2) (P 2 ) (P 2 ,A, (P 2 ( ^ <(-4 2 )r»> 1 _ 
P G^ip ,M ) = ^(2)^) - c ° (j^OO) + <* { 7 Mri) 4(jVg _ 1} ^ ■ (6) 

where G^{p 2 ) is the bare propagator. This expansion does not exactly separate the perturbative form the non- 
perturbative contribution because of the denominator /i 2 G^ (/j, 2 ) which contains a non perturbative contribution. It 
is therefore convenient to introduce a slightly different renormalisation R': 



p 2 G^(p 2 ) _ c o(^«(^)) <* (£,«*(/*)) <(A 2 ) R ( M )> i 



P K ' [P ' M } ~ v 2 G {2 l( P 2 ) co(l,a( M )) + c (f,a(M)) 4(JV? - 1) •* ' 



where (^4 2 ) R (fJ>) represents the A 2 operator renormalised as a local operator at the scale \i. Here the denominator is 
only the perturbative contribution to the Green function whence the first term in eq. (7) is purely perturbative: it 
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runs perturbatively with a pcrturbative MOM renormalisation condition at p 2 = p? . Let us define for simplicity the 
constant 

From [3] the first term in eq. (7), zo Co (p 2 / /J? , a;(/i)) , represents the three loop perturbative contribution and the sec- 
ond, z c 2 (p 2 / fi 2 , a(fj,)) , the leading logarithm Wilson coefficient of the 0(l/p 2 ) nonperturbative correction attributed 
to the vacuum expectation value of the renormalised local operator < (^4 2 ) R (/«) >• 
Let us now introduce some notations: 

(A 2 ) R ( M ) = Z- A l( P ) : A 2 arc := Z~\^ : (A R ( M )) 2 : , (9) 
where Z{p) = Z% (ii)Zjp{n), and 

™ < a <"» S dhV ln ^ (Al) = ^ + - ' 

-/ i xx _ rf 1 9/ x - . 3iVc oQa) 

7(a(M)) =dhv lnZ(M)= - 7o ^r + - = -- r~ + - (10) 



where the symbol : ... : represents the normal ordered product in the perturbative vacuum^. Our main goal in 



this paper is to understand better the connection between < (^4 2 ) R (m) > defined in [3] and the (A R (p)) 2 object 



considered here. 

The expansion in (7) is only valid above some momentum p > p m i n . Typically we have taken p m i n — 2.6 GeV for 
our fits reported in [1] - [4]. 

From cq. (5), (7) we decompose 



where 



< (A R (/i)) 2 > A =< (/i)) >r* + < {amy >r E + < (amy >i R (ii) 



< (A R (/i)) 2 >- rt = 3{N l~! )Z0 t dp 2 c ( 4,«(M)) ; (12) 



167T 2 J p 2 \ fj 2 

3(N 2 - l)zo f A2 dp 2 (p 2 \<(A 2 ) R ( M )> 



I A , ss* OPE 3(iV 2 -l)z r dp 2 (p 2 



A(N 2 1) 



(13) 

< (A R ( M )) 2 >? = ^-"il p 2 dp 2 G$ (p 2 ). (14) 

A few comments are in order here. < (A R (/x)) 2 >A Crt corresponds indeed to the perturbative computation of the 
vacuum expectation value of A 2 , i.e. to the connected diagrams with no external legs and with one A 2 inserted. 
However, the coupling constant and the gluon fields in the diagrams have been consistently renormalised at the scale 
fi. To leading order eq. (12) leads to 



< (Anirir > 



HN 2 - 1) 



A 



In 



TO 

A \ \ W 



Aqcd 




2 



Aqcd 

which diverges more than quadratically. Note that the dependence in p 2 ^ is subdominant. 



t^The : ... : symbols have been erroneously omitted in [3]. 
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In equation (13) the left hand side has been defined from the decomposition of the integral (5) according to (7). 
The right hand side contains < (^4 2 ) R (m) > already discussed. The latter is just a number which factorises out of 
the integral in (13). We thus see that < (t1r(^)) 2 >° PE an d < (^ 2 ) R (m) > are proportional. 

Our next task is to compute the proportionality coefficient and to compare < (Ar([i)) 2 >° PE with the other 
subleading term, < (A-^(p)) 2 >^. From [3] and eq. (8) we know that zq = 1 + 0(l//i 2 ). Our calculation of the 
integral in eq. (13) being performed to leading logarithm we will take zqC2 = C2 in the following. From (23) in [3], 



TO 



Let us also recall 



with 



TO /on T A 2 

Z 3 {p) cx (a( M ))^ , < (A 2 ) R ( M ) >cx (a^))"^ (17) 



35 9 
Po = 11, 7o = 13/2, 7,42 = — , 7o = 7A 2 - 7o = - (18) 



From eq. (13) and the leading logarithm relation 



47r da 

dp' lv" = rflog(p') ~ - 

we find 



d//p 2 =dlog(p 2 )~-^^ (19) 
Po " 



<m r,^ 2 . Q pE - g = 1} (127r) < (A2) ^ M > 4 " r (Pm,n) rf«r>- 1 - 

<( ^' (Ai)) >A -^6^—— i 4(iV? - 1) pb ' 

(a(/x)) ^ 



<(^ 2 ) R (M)> 



TO , TO 



a(A) \ ^ _ ( a(pmin) 



(20) 



It is interesting to notice that the coefficient pb/7o stemming from the integration over a is exactly compensated by 
the prefactors outside the integral, the origin of which does not appear at first sight to be related to the anomalous 
dimension of A 2 . Had we taken any other anomalous dimension instead of 70, say some 7', we would have ended with 
a constant 9/(47') in front of the r.h.s of eq. (20). 

In the large A limit, a(p m i n ) a(A) whence, since 70 is positive, the main result of this note comes from: 



TO 



«*•«•»' >*" a - 0o <WhW>(^)^- PU 

To leading logarithms and keeping \i fixed, 

< (Aji(fi)) 2 >° PE cx a(A)"7f A -> ^00 . (22) 

On the other hand, form eq. (14) 

< (A R (^)) 2 >^ r = constant. (23) 

since it does not depend on A. It results that < (Ar(p)) 2 >a PE is dominant over < (A-r([j,)) 2 >^ in the decomposition 
(11). This dominance will lead to the coming interpretation (next section) of the decomposition (11), which is indeed 
the main result of this note. Of course, we perform our analysis in the Landau gauge because of its conceptual [6] 
and numerical (Lattice Green functions [3]) particular interest. However the survival after renormalisation of BRST 
invariance in covariant gauges for the generalised composite operator showed in ref. [5] seems to point out that an 
analogous analysis, with similar results, for these gauges could be performed. 
As an interesting special case, if fj, = A 

< (A R (A)) 2 >° PE A - ^ < (A 2 ) R (A) >oc (a(A))"T*r (24) 
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III. CONCLUSION AND DISCUSSION 



Our conclusion is summarised in 



<(a%(»)> * < (a r ( M )) 2 >r E ( 



, 2fl 

a(A) - 



(M)) 



70 

< (A R (v)f > A - < (A R ^)) 2 >r] (^) *° • (25) 

The notations being not conventional let us recall that the < ... >a's in the r.h.s represent the gluon propagator 
integrated over momentum up to an UV cut-off, A, see eqs. (5) and (12). The gluon fields and coupling constants 
are renormalised in all the terms appearing in these equation. Thus we learn that the further renormalisation of the 
local operator A 2 proceeds by substracting to the plain vacuum expectation value of A 2 the same object computed in 
perturbation. This logarithmically divergent difference is then renormalised by the powers of a in the r.h.s of eq. (25). 
Not unexpectedly, we retrieve in essence the initial expression of the renormalisation of the A 2 operator through normal 
ordering (i.e. subtraction of the perturbative v.e.v), followed by the multiplicative, logarithmic renormalisation Z A 2. 
But apart from a non trivial consistency check, involving in particular the detailed expression of the Wilson coefficient, 
we obtain an expression which is more transparent, since it only involves a measurable quantity, the integral over the 
renormalised propagator. 

Eq. (25) presents a separation between perturbative and non-perturbative contributions to the integrated propaga- 
tor i.e. to A 2 . Of course, such a separation depends on the renormalisation scheme, and on the order in perturbation 
theory in which the Green functions are computed. It is also well known that summing to infinity the perturbative 
series may generate renormalons which behave like non-perturbative condensates. To avoid any such problem we 
stick to a finite order in the perturbative series. Furthermore, if the quantitative separation between perturbative 
and non-perturbative contributions depends on these prescription, the results summarised in eq (25) do not depend 
on them provided that we use the same scheme and order when computing both sides of eq. (25). Of course the 
anomalous dimensions to leading logarithms do not either depend on them. 

This simple result has several interesting consequences. First, it has been advocated [4] that the A 2 condensate 
could be dominantly due to the contribution to the path integral of semi-classical gauge field configurations such as 
instantons liquids. It is useful to consider this hypothesis through a background field picture, i.e. factorising the path 
integral into an integral over semi-classical gauge field configurations, and for each value of these an integral over 
quantum fluctuations around this background configuration. It means that the hermitian matrix A^ is decomposed 
into: 

A^ = B^ + Q^B) A 2 =B 2 + {B.Q}+ + Q 2 (B) (26) 

B^ being the background, assumed to be non-perturbative, and the quantum fluctuations assumed to be pertur- 
bative. {B.Q} + = B.Q + Q.B. In principle, depends on B^ and differs from the quantum fluctuations around 
the trivial vacuum B^ — which is what perturbative QCD computes. The hypothesis that < (^4 2 ) R > is due ^ to 
these semi-classical gauge configurations is translated into: (^4 2 ) R — B 2 . From eqs. (25), (26) 

< (A 2 ) R >~< B 2 >~< A 2 > - <Q 2 (B = 0) > (27) 

i.e. that (Q 2 (B) — Q 2 (B = 0)) is subleading §§ . In other words the dependence of Q M on B^ is subleading. The hard 
quantum fluctuations are not sensitive to the soft background field. 

A most interesting consequence of our result is related to some discussions in [6]. These authors extend to QCD 
some remarks stemming from compact £7(1). They attribute a special role to the A 2 condensate, even if a gauge 
dependent quantity, by arguing that A 2 in the Landau gauge is the minimum of A 2 on the gauge orbit. One difficulty 



"This discussion is qualitative and we do not know how to define rigorously the corresponding scale /i. We therefore prefer 
to omit writing /i here. 

§§ If B is a classical solution of the field equations, the term linear in Q will vanish. B should be close to such a solution and 
we therefore neglect {B.Q} + . 
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in this argument is the following: Fixing the Landau gauge amounts to minimize the < > while the condensate 

refers to some renormalised quantity free of the quadratic and logarithmic divergences. In compact U(l) life is simpler: 

< ^barc > = < ^pcrt > + < ^nonpert >' (28) 

the pcrturbative theory is trivial and the nonperturbative contribution is due grossly speaking to the topology. A phase 
transition when the coupling constant varies allows to measure directly the non perturbative contribution. We refer 
to [6] for more details. Our result eq. (25) exhibits in QCD, up to subleading contributions, a linear decomposition 
similar to eq. (28), although such a similarity is not at all obvious at first sight. The next question could be whether 
in some sense the < A 2 > OPE computed in the Landau gauge is the minimum of some quantity on the gauge orbit. 

Last but not least, let us simply say that the result in eq. (25) provides a fairly simple understanding of what the 
A 2 condensate is. It confirms that indeed the 0(\/p 2 ) correction to perturbative QCD at large momenta has to do 
with the A 2 condensate. Indeed, if one start with some doubt about the relation of the r.h.s of eq. (13) with an A 2 
condensate, just considering it as an unidentified 1/p 2 contribution, we end-up with the conclusion that it yields a 
non-perturbative contribution to the A 2 v.e.v. The fact that in our derivation this term has precisely the anomalous 
dimension of an A 2 condensate comes form the fact that c 2 in the r.h.s of eq. (13) has been computed under the 
assumption that it is due to an A 2 condensate, an assumption which has been shown to fit fairly well the lattice data. 
Had we used another scale dependence for c 2 we would have ended with a wrong scale dependence for the resulting 
non-perturbative contribution to the A 2 v.e.v. We would have also ended with a constant different from 1 in front of 
the r.h.s of eq. (25), see the discussion following eq. (20). Thus the picture is fully consistent. 
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